J1128+592: a highly variable IDV source 
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The dates of receipt and acceptance should be inserted later 
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Short time-scale radio variations of compact extragalactic radio quasars and blazars known as IntraDay Variability (IDV) 
can be explained in at least some sources as a propagation effect; the variations are interpreted as scintillation of radio 
waves in the turbulent interstellar medium of the Milky Way. One of the most convincing observational arguments in 
favor of a propagation-induced variability scenario is the observed annual modulation in the characteristic time scale of 
the variation due to the Earth's orbital motion. So far there are only two sources known with a well-constrained seasonal 
cycle. Annual modulation has been proposed for a few other less well-documented objects. However, for some other IDV 
sources source-intrinsic structural variations which cause drastic changes in the variability time scale were also suggested. 
11128+592 is a recently discovered, highly variable IDV source. Previous, densely time-sampled flux-density measure- 
ments with the Effelsberg 100-m radio telescope (Germany) and the Urumqi 25-m radio telescope (China), strongly in- 
dicate an annual modulation of the time scale. The most recent 4 observations in 2006/7, however, do not fit well to the 
annual modulation model proposed before. In this paper, we investigate a possible explanation of this discrepancy. 
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1 Introduction ing arguments in favor of an extrinsic explanation of IDV 

is the so called annual modulation of the IDV time scale. 
Short time-scale (few hours to few days long) variations of The characteristic variability time scale is inversely propor- 
the radio flux density in flat spectrum radio-l oud quasars t ional to the relative velocity between the observer and the 
and b l azars were discove red in the mid-eighties dHeeschen et alt ering medium. The observer's velocity (and so the rel- 
1987HWitzeletalJll986l) . The phenomenon was named In- ative velocity vector between the observer and scattering 



medium) undergoes a systematic annual modulation as the 
Earth orbits around the Sun. This annual velocity variation 
then is observed as an annual change in the variability time 
scale. Such s easonal cycles are seen in two ID V sources: 
J181 9+3845 jPennett-Thorpe & de Bruvnl2003l) and PKS1257-| 
326 (iBignall et al. l2003|). In a few oth er IDV sources, such 



traDay Variability (IDV). If interpreted as being source in- 
trinsic, the short time scale of the variations would imply - 
through the light travel time argument - micro-arcsecond- 
scale sizes of the emitting regions, which would result in 
excessively large apparent brightness temperatures. The cal- 
culated variability brightness temperatures obtained are typ- 
ically in the range of 10 16 K to 10 21 K, wh ich is far in excess a s PKS 15 19-273 djauncev et al.ll2003l). 
of the inverse-Compton limit of ~ 10 12 K dKellermann & Pa ulirff 0fdtn|624 djauncev & Macauard2001r.lRickett et all200lhl 

PKS0405-385 dKedziora-Chudczerll2006l) and B0954+658 
(Fuhrmann priv. com.), different characteristic variability 
time scales were measured at different observing epochs, 
however systematic seasonal variations, as expected from 
the annual modulation model, were either not seen or could 



1969). Thus, theories which explain IDV with variations 
intrinsic to the quasar, require excessively large Doppler 
boosting factors (6 » 50) or special source geometries (such 

as n on-spherical rela t ivistic emission mode ls, 

e.g. lOian et al]fl99ll Il996albl ISoada etal.l[l 999) or coher- 



ent an d collective plasma emission dBenfordll992r.lLesch & PoKlf be unam biguously identified. T he IDV sources 



1992) to avoid the inverse-Compton catastrophe. 

An alternative theory explains IDV as a propagation ef- 
fect. In this source-extrinsic interpretation, IDV is caused 
by interstellar scintillation (ISS) of radio waves in the tur- 
bulent plasma of the Milky Way. One of the most convinc- 



B0917+624 dFuhrmann et alJ2002l) and PKS0405-385 seem| 
to display the so called e pisodic IDV (e.g. 
Kedziora-Chudczerll2006l) . where pronounced variability is 
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observed for months or even years and then temporarily 
stops and the IDV ceases. This variability behaviour may be 
due to changes in the intrinsic source structure (expansion 
of the scintillating component) or changes in the scatter- 



'*WILEY 

jinterScience* 



© 2006 WILEY- VCH Verlag GmbH & Co. KGaA, Wcinheim 



ing plasma (Bernhart et al . 2006; Ked ziora- Chudczerl2 006: 



Kraus et al.lfl999tlKrichbaum et alj|2002l 

If IDV is caused by ISS, it is possible to predict the 
frequency dependence of the variability amplitude, which 
above a transition frequency from strong to weak scintilla- 
tion is expected to be less pronounced towards higher ob- 
serving frequencies. However, obs ervations of several IDV 
sources (most notably S50716+71. lAgudo et alJ2 006: Kraus 
12003b lOstorero et alJ 120061) do not show such a frequency 
dependence and show variability amplit udes which increase 
with frequency. iKrichbaum et all (2002) suggested that a 
source-intrinsic contribution may cause the discrepancy be- 
tween the observed and the expected frequency dependence. 

For a better understanding of the IDV phenomenon, a 
clear distinction between ISS induced propagation effects 
and a source-intrinsic origin of the variability is important. 
Here, we present the variability characteristics of the re- 
cently discovered highly variable IDV source J 1128+592 
and discuss the possible mixture of ISS-induced and source- 
intrinsic variability. 

2 Observations 

J 1128+592 was observed with the Max-Planck-Institut filr 
Radioastronomie (MPIfR) 100-m Effelsberg radio telescope 
(at 2.70, 4.85 and 10.45 GHz) and with the Urumqi radio 
telescope (at 4.85 GHz) in so far 15 observing sessions, each 
lasting several days. The results of the first ten epochs o f 
observations are already published in lGabanvi et a l] (l2007l) . 

The low-noise 4.85 GHz receiver, a new receiver back- 
end and new telescope driving software for the Urumqi tele- 
scope were provided by the MPIfR. A detailed technical de- 
script ion of the recei ver system and the telescope is given in 
e.g. lSunetai1(l2007b . 

The observing epochs are irregularly distributed over 
a ~ 2-year interval between December 2004 and January 
2007. The 4.85 GHz observing epochs, and the observing 
radio telescopes are listed in Col. 1 and Col. 2 of Table [T] 

At both telescopes, the flux density measurements were 
performed using the cross-scan technique, in which the tele- 
scope is moved repeatedly (4-8 times) in azimuth and 
in elevation over the source position. Each such movement 
is called a subscan. After baseline subtraction, a Gaussian 
curve was fitted for each subscan to the resulting slice across 
the source. For each scan, the averaged and pointing cor- 
rected peak amplitude of the Gaussian curve yielded a mea- 
sure of the source flux density. Then the systematic eleva- 
tion and time-dependent gain variations were corrected, us- 
ing the combined gain curves and gain-transfer functions 
of calibrator sources of known constant flux density. The 
measured flux densities were then tied to the absolute flux- 
density scale. The absolute flux-density scale was deter- 
mined from repeated observations of the primary calibra- 
tors e.g. 3C48, 3C286, 3C295 and NGC 7 027 and using 
the flu x-density scale of Baars et al. ( 1977 ) and Ott et al.1 
dl994l) . 



Table 1 Summary of the 4.85 GHz IDV observations of 
Jl 128+592. The table lists the observing dates (Col. 1), the 
observing radio telescopes (Col. 2, 'E' for Effelsberg, 'U' 
for Urumqi), the average modulation index of the calibrators 
(too, Col. 3), the modulation index of Jl 128+592 (Col. 4) 
and the derived characteristic variability time scale (Col. 5) 
in days. 
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a iGabanvi et all 120071) 

'' There were not enough measurements to calculate a reliable time 
scale. 

c Due to human error, slightly different frequency setup was 
used than previously. The center frequency was 4.79 GHz instead of 
4.85 GHz. 



A more detailed description of the data reduction and 
the compa rison of accuracy rea ched by the two telescopes 
is given in lGabanvi et alJ ( 2007 ). 



2.1 Data analysis 

For the variability analy s is, we follow the method intro- 
duc ed by Heeschen et alJ ( 1987[) and desc r ibed i n detail in 
e.g. lOuirrenbach et all d2000l) , lKraus et alJ d2003l) . 

To decide whether a source can be regarded as variable, 
we performed a x 1 test- We tested the hypothesis whether 
the light curve can be modeled by a constant function. The 
sources for which the probability of this hypothesis was less 
then 0.1% were considered to be variable. The amount of 
variability in a light curve is described by the modulation in- 
dex and the noise-bias corrected variability amplitude. The 
modulation index is defined as to = 1 00 ■ ^ , where (S ) is the 
time average of the measured flux density, and <x is the stan- 
dard deviation of the flux density. The average modulation 
index of the calibrator sources (too) provides a conservative 
estimate of the overall calibration accuracy. It describes the 
amount of residual scatter in the data set after all the cali- 
bration steps are performed. The too and to of J 1128+592 
are listed in Col. 3 and in Col. 4 of TableQ] respectively. 
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To determine the characteristic variability time scale, 
we m ade use of the struc ture function (for a definition see 
e.g. ISimonetti et alJI 19851) . the autocorrelation function and 



the light curve itself. We defined the characteristic variabil- 
ity time scale by the time-lag, where the structure func- 
tion reaches its saturation level. This corresponds to the first 
minimum of the autocorrelation function. We compared 
these with the average peak- t o- trough time derive d from the 
light curve (for details see IGabanyi et al.l 120071) . In Col. 5 
of TableQ] we list the variability time scales obtained by the 
average from the three different methods. 

3 Results and discussion 

In Fig.[T|the light curves of five new, previously unpublished 
IDV observations of J 1128+592 are shown. As in all the 
previous measurements, J 1 128+592 showed pronounced vari- 
ability with peak-to-trough amplitudes of up to 20-25 %. 

The observed variability time scales of J 1128+592 are 
signif icantly differen t at dif ferent dates (see Col. 5 of Table 
[T]). In iGabanyi et al.l d2007l) . we proposed that the changes 
in the characteristic variability time scale may be due to an- 
nual modulation. The best fit to the data was a chieved using 
the an isotropic annual modulation model of Bignall et alJ 
d2006l) . which is represented by the solid curve in Fig. [2] 

In the anisotropic model, the scintillation time scale also 
depends on the ellipticity of the scintillation pattern and on 
the direction in which the relative velocity vector (between 
the Earth and the screen) "cuts through" the elliptical scin- 
tillation pattern. Thus, the fitted parameters obtained from 
the anisotropic scintillation model are the velocity compo- 
nents of the scattering screen, the scattering length-scale 
(which depends on the screen distance and the scattering 
angle), the angular ratio of the anisotropy and its position 
angle. 

However, not all of the new observations (data points at 
day 234, day 267, day 309 and day 13) fit equally well to this 
model. It is clear that the time scales derived from the ob- 
servation in September 2005 and in September 2006 do not 
agree well with each other, questioning the interpretation 
via annual modulation. The larger error bar on the time scale 
obtained in September 2006 is due to the fact that only one 
well-defined variability peak occurred during this observa- 
tion. Moreover the relatively slow time scale as measured 
in November 2006 (day 309) does not agree well with the 
model. On the other hand, the acceleration of the variability 
time scale, which is expected after autumn (day > 240) is 
basically seen again, also with the new data. Despite some 
difficulties in the determination of the time scale in Decem- 
ber 2006 (due to irregular sampling), the variability in this 
months is faster than in autumn. Also the measurements of 
January 2007 (day 13) are consistent with such fast vari- 
ability, although this recent data point shows a bit slower 
variations than in February 2006 (see Fig. [TJ. 

Inspecting the structure function and autocorrelation 
function of the 2006 November data reveals a possible faster 



time scale of ~ 0.5 day. This time scale would fit to the 
model. Unfortunately the sampling and the amount of mea- 
surements obtained in the later epoch, in 2006 December 
are not adequate to check the existence of two time scales. 
At the previous epoch (2006 September), there is no sign of 
variation on an additional time scale either. More observa- 
tions are necessary to confirm or reject the appearance of a 
secondary time scale. 

Apart from the time scale, the modulation index shows 
significant variations as well. Notably, in 2006 September 
the modulation index is half of the one measured in the pre- 
vious epoch (2006 August) and the following epoch (2006 
November). In the latter case m = 4.7 %, whilst in Septem- 
ber m = 2.3 %. The annual modulation model cannot ex- 
plain changes in the strength of variability. These different 
variability indices might therefore suggest changes in the 
-scattering plasma or intrinsic changes in the source. If the 
scattering angle or the intrinsic source size has increased 
in this epoch, we would expect a prolongation of the time 
scale as well. So this might explain the discrepancy with the 
annual modulation model. 

Finally, we note that the mean flux density at 4.85 GHz 
changed significantly during the ~ 2 years of monitoring 
of Jl 128+592. The flux density increased by ~ 27 % until 
2006 February, then monotonously decreased until our last 
observation to approximately the same flux-density level. 
This long term flux density change suggests a source-intrin- 
sic origin (e.g. an ejection of a new jet component), which 
could influence the variability behavior of the source on 
IDV time scales as well. Already proposed Effelsberg mon- 
itoring observations will help us to further investigate these 
questions. Additionally, proposed Very Long Baseline Ar- 
ray observations can reveal any intrinsic changes in the ra- 
dio structure of Jl 128+592. 
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Figure 1 4.85 GHz light curves of Jl 128+5925 (upper plots) and those of the secondary calibrator B0836+710 (lower 
plots) from 2006 August (panel (a) ) until 2007 January (panel (e) ). On the abscissa the observing time (in days) is 
displayed, the starting dates are given in the upper right corner of each plot. On the ordinate are the normalized flux 
densities in Jansky. Panels (a), (b) and (c) show Urumqi observations, panels (d) and (e) show Effelsberg observations. 
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Figure 2 The characteristic variability time scales of 
Jl 128+592 plotted versus day of the year. The data are 
from 4.85 GHz observations performed with the Effelsberg 
telescope (filled symbols) and with the Urumqi telescope 
(open symbols). Different symbols represent observations 
performed in different years: star stands for 2004, squares 
for 2005, circles for 2006 and triangle for 2007. The solid 
curve re presents the best fit annual modulation model as 



given in lGabanyi et a l. (2007), which uses only the first ten 



observing epochs for the fit. 
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